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Abstract 

The aim of this work is to unravel a basic microscopic picture behind complex magnetic prop- 
erties of hexagonal manganites. For these purposes, we consider two characteristic compounds: 
YMnO,3 and LuMnC>3, which form different magnetic structures in the ground state (P6 ? cm and 
PG^cm, respectively). First, we establish an electronic low-energy model, which describes the be- 
havior of the Mn 3d bands of YMnC>3 and LuMnOs, and derive parameters of this model from the 
first-principles calculations. From the solution of this model, we conclude that, despite strong frus- 
tration effects in the hexagonal lattice, the relativistic spin-orbit interactions lift the degeneracy of 
the magnetic ground state so that the experimentally observed magnetic structures are successfully 
reproduced by the low-energy model. Then, we analyze this result in terms of interatomic mag- 
netic interactions, which were computed using different approximations (starting from the model 
Hamiltonian as well as directly from the first-principles electronic structure calculations in the 
local-spin-density approximation) . We argue that the main reason why YMn03 and LuMnOs tend 
to form different magnetic structures is related to the behavior of the single-ion anisotropy, which 
reflects the directional dependence of the lattice distortion: namely, the expansion and contraction 
of the Mn-trimers, which take place in YMnC>3 and LuMnC>3, respectively. On the other hand, 
the magnetic coupling between the planes is controlled by the next-nearest-neighbor interactions, 
which are less sensitive to the direction of the trimerization. In the P63 cm structure of YMnC>3, 
the Dzyaloshinskii-Moriya interactions lead to the spin canting out of the hexagonal plane - in the 
same direction as the single-ion anisotropy. Finally, using the Berry-phase formalism, we evaluate 
the magnetic-state dependence of the ferroelectric polarization, and discuss potential applications 
of the latter in magnetoelectric switching phenomena. 
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I. INTRODUCTION 



Hexagonal manganites (the space group PQ^cm) are one of canonical examples of multi- 
ferroics, which have attracted an enormous attention recently. The coexistence of ferroelec- 
tricity and magnetism in such systems provides a unique possibility for manipulating the 
charges by applying a magnetic field and the spins by applying a voltage, which is crucially 
important for the construction of new forms of multifunctional devices.- To this end, the 
direct magnetic phase control by static electric field was realized in HoMn03~ The interplay 
between the ferroelectric activity and the magnetic order was also demonstrated in YMnC>3 
and LuMnOs with the measurements of the dielectric constant and the loss tangent, which 
were shown to exhibit clear anomalies around the Neel temperature (Tn= 75 K and 88 K in 
YMnOs and LuMn03, respectively),- 1 ^ even despite the fact that the ferroelectric transition 
itself occurred at much higher temperature (Tc~ 880 K). 5 Another spectacular example is 
the coupling of magnetic and ferroelectric domains, which was visualized in YMnC>3 by us- 
ing optical second harmonic generation technique.- Furthermore, the magnetic transition in 
YMnC>3 and LuMn03 is accompanied by a distinct change of the atomic positions. 7 - Thus, the 
experimental data clearly demonstrates the existence of a strong coupling amongst electric, 
magnetic, and lattice degrees of freedom in these hexagonal manganites. 

The magnetic frustration is one of the key concepts of multiferroic materials, which may 
assist the inversion symmetry breaking and, in a number of cases, be even responsible for 
such a breaking. 8 In this respect, the hexagonal lattice is not an exception, and is typi- 
cally regarded as a playground for studying the magnetic frustration effects. However, it is 
also the main complication, hampering the theoretical understanding of multiferroic effects 
in hexagonal compounds, even despite the fact that the high-spin state (S=2), realized in 
manganites, is typically regarded as an "easy case" for such theoretical analysis, where the 
classical the spin fluctuations dominate over the quantum ones. Nevertheless, the ground 
state of classical spins in the hexagonal lattice is expected to be highly degenerate. Different 
signs of spin fluctuations, apparently originating from this degeneracy, were indeed observed 
in the neutron scattering experiments, even below T^.—^ Another evidence of spin fluctua- 
tions, which is also related to the quasi-two-dimensional character of magnetic interactions, 
is the large ratio of the Curie-Weiss temperature (#cw) to Tn (about 7 in YMnOs).- 

The degeneracy can be lifted by lattice distortions and, in this context, plenty of attention 
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is paid to the so-called trimerization instability, inherent to the P'o^cm structure.-^ How- 
ever, the trimerization alone does not lift the frustration of isotropic exchange interactions. 
In this sense, the situation is fundamentally different from the exchange striction effect, 
which accompanies the formation of the E-type antiferromagnetic (AFM) state in the or- 
thorhombic YMnC>3 and which lifts the frustration of nearest-neighbor (NN) interactions.— 
Nevertheless, the trimerization can interplay with the relativistic spin-orbit (SO) coupling 
and, in this way, give rise to new anisotropic interactions, which can lift the degeneracy and 
stabilize some individual magnetic structure with the well-defined symmetry Such struc- 
tures were detected in the experiments on the neutron diffraction (Ref. 
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M) and 



optical second harmonic generation (Ref. Il5l). In a number of cases (e.g., in LuMnOs), there 
can be several magnetic structures, coexisting in a narrow temperature range.— In short, 
despite difficulties, there is an enormous experimental progress in the identification of mag- 
netic structures of hexagonal manganites, resulting from a delicate balance between lattice 
distortion, SO interaction, and frustration effects. 

The microscopic understanding of rich magnetic properties of the hexagonal manganites 
is still rather limited. To begin with, there is no clear microscopic model, which would 
explain the origin of basic magnetic structures of hexagonal manganites, and why different 
manganites tend to form different magnetic structures. Basically, it is only known how the 
trimerization affects the NN isotropic interactions.— The presence of single-ion anisotropy 
and Dzyaloshinskii-Moriya (DM) interactions is, of course, anticipated. However, it is abso- 
lutely not clear how all these effects come together to form a variety of magnetic structures, 
realized in the hexagonal manganites. 

In this paper, we will try to answer some of these questions. For these purposes, we con- 
sider two characteristic manganites: YM11O3 and LuMnOa, which form different magnetic 
structures in the ground state: PQ^cm and PQ^crn, respectively (in the International nota- 
tions, where each underlined symbol means that given symmetry operation is combined with 
the time inversion). We will show that this difference can be naturally related to different 
directions of the trimerization: expansion and contraction of the Mn-trimer, which takes 
place in YM11O3 and LuMnOs, respectively. In our study, we start from the first-principles 
electronic structure calculations. First, we construct a low-energy electronic model, which 
captures details of the magnetic structure and correctly reproduces the magnetic ground 
state of YMn03 and LuMnOs. Then, we analyze these results by further transforming the 



electronic model into the spin one and elucidating which magnetic interaction is responsible 
for each detail of the magnetic structure. We will also consider the 'temperature effect', as- 
sociated with the temperature change of the experimental crystal structure, and show that 
above Tn it gradually diminishes the anisotropic interactions. 

The rest of the paper is organized as follows. All methodological aspects, such as con- 
struction of the electronic model and calculation of magnetic interactions, are discussed in 
Sec. [TTJ Results of solution of the electronic model in the Hartree-Fock (HF) approximation 
are presented in Sec. IIII Al In Sec. IIII B\ we give a detailed analysis of the obtained results 
in terms of magnetic interactions, which were computed using different starting points. In 
Sec. IIII C\ we discuss the magnetic part of the ferroelectric polarization and propose how it 
can be controlled by switching the magnetic state. Finally, a brief summary of the work is 
given in Sec. [IVl 



II. METHOD 



Since our goal is the construction of microscopic theory for the magnetic properties of 
YMnC>3 and LuMnOs, we first adopt the low-energy model, which would provide a realistic 
description for the Mn 3c? bands of these compounds: 

ij a/3 i {a} 

The model is constructed in the basis of Wannier orbitals, using the input from the first- 
principles electronic structure calculations. Each Wannier orbital is denoted by the Greek 
symbol, which itself is the combination of spin (s= f or \) and orbital (m= xy, yz, 3z 2 —r 2 , 
zx, or x 2 —y 2 ) variables. Since the Mn 3c? bands in hexagonal manganites are well separated 
from the rest of the spectrum,— the construction of the model Hamiltonian ([!]) is rather 
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straightforward. The corresponding procedure can be found in Ref. 

All calculations have been performed using experimental parameters of the crystal struc- 
ture, measured at 10 K and 300 K (Ref. 12, Supplementary Information), i.e. well below and 
above the magnetic transition point. The experimental space group PQ^cm has 12 symmetry 
operations, which can be generated by the mirror reflection x— Y— x, m x , and the 60°-degree 
rotation around the z-axis, combined with the half of the hexagonal translation, {C||c/2}. 

The crystal-field splitting, obtained from the diagonalization of the site-diagonal part 



of Uj = • || (without spin-orbit coupling), is very similar in YM11O3 and LuMnOs. For 
example, if one uses the 10 K structure, we obtain the following scheme of the atomic levels: 
-0.54, -0.43, -0.29, -0.24, and 1.50 eV in the case of YMn0 3 , and -0.60, -0.49, -0.25, 
—0.24, and 1.57 eV in the case of LuMnOs. The use of the 300 K structure yields similar 
results. Clearly, the crystal field tends to stabilize four atomic orbitals, which are separated 
from the fifth one by the large energy gap. Such a scheme of the crystal-field splitting is 
consistent with the formal d A configuration of the Mn-ions, which is subjected to the Jahn- 
Teller instability. The fifth (unoccupied) orbital is of predominantly 3z 2 —r 2 symmetry. The 
off-diagonal elements of iij = \\t"- | w ith respect to the site indices stand for the transfer 
integrals. They are listed in Ref. 171 The value of the screened Coulomb interaction U 
(defined as radial Slater's integral F°) is about 2.6 eV for all considered systems. The 
intraatomic exchange (Hund's) coupling Jh is about 0.9 eV, which is practically unscreened. 
The full matrices of screened Coulomb interactions can be also found in Ref. [17. 

After the construction, the model is solved in the HF approximation.— This procedure 
appears to be extremely useful, especially for the search of the magnetic ground state. 
Typically, in frustrated magnetic systems, we are dealing with the competition of several 
magnetic interactions of the both relativistic and non-relativistic origin. Therefore, even HF 
calculations for the relatively simple model ([1]) can be very time consuming, because they 
may require tends of thousands of iterations. In such a situation, the full scale electronic 
structure calculations are simply unaffordable. Since the degeneracy of the ground state is 
lifted by the lattice distortion, the HF approximation appears to be a good starting point 
for the analysis of the equilibrium magnetic properties.— 

Of course, the model (CEJ) is not perfect, because it does not explicitly include the oxy- 
gen band, which can be important for the quantitative analysis of magnetic properties of 
manganites. Therefore, whenever possible, we check results of our model analysis by di- 
rect calculations in the local-spin-density approximation (LSDA). For these purposes, we 
use the tight-binding linear muffin-tin-orbital method (in the following we will refer to such 
calculations as 'LMTO calculations').— Hopefully, in both cases we can employ the same 
strategy for calculations of magnetic interactions, which is based on the local force theorem 
and the Green's function technique. Namely, the isotropic exchange interactions (Jij) can 
be obtained in the second order perturbation-theory expansion for the infinitesimal spin 
rotations,— antisymmetric DM interactions (<%) - by considering mixed type perturba- 



tion with respect to the rotations and the relativistic SO coupling,—"— and the single-ion 
anisotropy tensors (fj) - in the second order with respect to the SO interaction^. 

The LMTO calculations have been performed for the AFM configuration tititi) where 
the arrows stand for the directions of magnetic moments at the sites 1-6 (see Fig. [TJ for 
the notations of atomic positions). The use of the AFM configuration is essential in order 
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FIG. 1. (Color online) Relative positions of Mn-sites in the hexagonal PQ^cm structure: the atoms 
located in the plane z=0 are indicated by the red (dark) spheres, and the atoms located in the 
plane z=c/2 are indicated by the light orange (grey) spheres. The Mn-trimers, which transform to 
each other by the symmetry operation {Cf |c/2}, are shaded. 

to open the band gap in LSDA (about 0.7 eV for YMnOa, which is comparable with the 



experimental optical gap of 1.3 eV, reported in Ref. 1241 ) . Certain inconvenience of working 
with the AFM tititi configuration is that it artificially lowers the PQ^cm symmetry: in 
this case, the local symmetry can by preserved only around the sites 2 and 5, which will be 
selected as the reference points for the analysis of interatomic magnetic interactions. 

In our LMTO calculations we decided to stick to the regular LSDA functional and not to 
use any corrections for the on-site Coulomb interactions (LSDA+[7). On the one hand, such 
corrections can improve the description for interatomic magnetic interactions (similar to the 
model). On the other hand, the use of the LSDA+?7 functional is always conjugated with 
some additional uncertainties in the calculations, related to the double- counting problem. 



Furthermore, the example of LaMn0 3 shows that LSDA is a reasonably good starting point 
for the analysis of interatomic magnetic interactions.— Nevertheless, when we compare the 
LMTO results with the model calculations we discuss possible consequences of the Coulomb 
U on the magnetic interactions in the former case. 

Due to the hybridization with the oxygen states, which is treated explicitly in the LMTO 
calculations, the value of spin magnetic moment at the Mn-sites is reduced till 3.5 /xb- 
Thus, some deviation of the local magnetic moment from the ionic value (4 /ib), which is 
typically seen in the experiment,— 1 ^ can be attributed to the covalent mixing. In model HF 
calculations, similar effect can be described through the transformation from the Wannier 
basis to that of atomic orbitals. 16 



III. RESULTS AND DISCUSSIONS 



A. Optimization of Magnetic Structure 



We start with the central result of our work and argue that the low-energy model ([I]), 
with the parameters derived from the first-principles electronic structure calculations,— suc- 
cessfully reproduces the magnetic ground state of YMn0 3 and LuMn0 3 . 

The main candidates for the magnetic ground state of YMn0 3 and LuMn03 are shown 



in Fig. |2] (see also Refs. |l3| and [l4j for the notations). The unidimensional representations 
Ti, T2, Vs, and T4 correspond to the magnetic space groups P6 3 cm, P6 3 cm, PQacm, and 
P6 3 cm, respectively. The directions of the spin magnetic moment, obtained in the HF 
calculations for the low-energy model, are listed in Table [H In the T\ and T4 configurations, 
all magnetic moments lie in the x?/-planes, while in the T2 and r 3 configurations, there is 
also a small canting of magnetic moments along z. Moreover, the T 2 configuration allows 
for the weak ferromagnetism along z, while in the T 3 configuration, the ^-components of 
the magnetic moments in the planes z=0 and z=c/2 cancel each other. More generally, the 
configurations Ti (T2) and T4 (r 3 ) differ by the magnetic alignment in adjacent xy-planes: 
{Cf|c/2} acts as the normal symmetry operation in Ti and which transforms these 
states to themselves, while in r 3 and T4, {C!?|c/2} enters the magnetic symmetry group in 
the combination with the time-inversion operation T. It corresponds to the additional flip 
of magnetic moments in the odd xy-planes of T 3 and T 4 . We have also considered other 
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FIG. 2. (Color online) Magnetic structures obtained in the calculations (in the notations of Ref. 

ri (a), r 2 (b), r 3 (c), r 4 (d), r 5 with ei||[ioo] (e), r 5 with ei ||[i20] (f), r 6 with ei ||[ioo] ( g ), 

and Tq with ei||[120] (h). The oxygen atoms are indicated by the small green (grey) spheres. The 
manganese atoms are indicated by the big spheres: the ones located in the z=0 plane are shown 
by the red (dark) color, and the ones in the z=c/2 plane - by the light orange (grey) color. 



possible magnetic configurations with the symmetries T 5 and r 6 , as explained in Ref. Il3 
However, as it will become clear from the discussion below, they have higher energies. 

The total energies of different magnetic configurations are summarized in Table [III 
Thus, the ground state of YMnC^ is T 3 (PQ 3 cm), in agreement with the experiment.— 1 ^ In 
LuMnOs, the ground state changes to I\ {PQ 3 cm), also in agreement with the experiment.— 1 ^ 
However, all the states are located in a narrow energy range, which is expected for frus- 
trated magnetic systems. The lower-symmetry magnetic structure PQ 3 , which is typically 
regarded as another possible candidates for the magnetic ground state of the hexagonal 
manganites,— appears to be unstable and steadily converges to either PQ 3 cm (YMnC^) 
or PG^cm (LuMnOs). 

The band gap, obtained for both YMnOs and LuMnOs, is about 2 eV, which is larger 
than experimental 1.3 eV.— Nevertheless, such a discrepancy is quite expectable for the level 
of the HF calculations. 
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TABLE I. The angles a and 0, representing the directions ej = (cos a« cos /3j, cos ai sin/3j, sinc^) 
of the spin magnetic moments in the plane z=0, for different magnetic configurations (results of 
calculations using the experimental crystal structure, measured at 10 K). The atomic positions are 
explained in Fig. [TJ For the magnetic configurations T%, T 2 , and Tg, the directions of the magnetic 
moments at the sites 4, 5, and 6 in the plane z=c/2 are obtained by the 180° rotations around the 
z-axis of the vectors ei, 62, and e3. For the magnetic configurations T%, T4, and T§, these 180° 
rotations should be combined with the time inversion. 
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B. Analysis of Magnetic Interactions 

In this section, we clarify results of the HF calculations for the low-energy model and argue 
that such a good agreement with the experimental data for the magnetic ground state is not 
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TABLE II. Total energies of different magnetic configurations as obtained in the Hartree-Fock 
calculations for the low-energy model. The energies are measured in meV per one formula unit, 
relative to the most stable configuration. The magnetic configurations are explained in Fig. [2J 
The calculations for YMn03 and LuMn0 3 have been performed using the experimental crystal 
structure, measured at 10 K and 300 K (as denoted in the notations). 

configuration YMn0 3 (10 K) YMn0 3 (300 K) LuMn0 3 (10 K) LuMn0 3 (300 K) 
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surprising and can be understood from the analysis of corresponding magnetic interactions, 
which in turn depend on details of the lattice distortions in YMn0 3 and LuMn0 3 . Thus, 
we consider the spin model: 

Hs = - ^2 ■ i >j°- ;v 'j + di ^ e{ x e ji + e ^ e *> ( 2 ) 

(ij) (ij) i 

which can be obtained by eliminating the electronic degrees of freedom form the more general 
Hubbard model (00), or directly from the LMTO calculations.—^— In these notations, {Jij} 
are the isotropic exchange interactions, {<%} are the antisymmetric DM interactions, {fj} 
are the single-ion anisotropy tensors, stands the direction of the spin magnetic moment 
at the site i, and the summation runs over all pairs of atoms (ij). 

The parameters of isotropic magnetic interactions are listed in Table IHIt and the atomic 
positions are explained in Fig. [T] All NN interactions in the plane xy are AFM. This 
is reasonable, because the ferromagnetic (FM) coupling in the hexagonal geometry can be 
stabilized only by virtual hoppings onto unoccupied 3z 2 —r 2 orbital, which are relatively small 



(see Ref. 



171 ). Moreover, the number of orbital paths, available for the virtual hoppings via 



this particular 3z 2 — r 2 orbital, is also small. Nevertheless, from the orbital decomposition of 
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TABLE III. Parameters of isotropic exchange interactions (measured in meV), calculated in the 
ferromagnetic states of YMn0 3 and LuMn03. The atomic positions are explained in Fig. [TJ 
Calculations have been performed using the experimental parameters of the crystal structure, 
measured at 10 K and 300 K (as denoted in the notations). 
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-0.20 
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-0.19 


-0.11 


-0.08 
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-0.24 


-0.31 


-0.21 


-0.24 
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-0.07 


-0.16 


-0.16 


-0.23 



Jij in our LMTO calculations, we can conclude that such a FM contribution does exists and 
compensates about 30 % of AFM contributions, involving all other orbitals, except 3z 2 —r 2 . 

The symmetry of the PQ^cm lattice is such that there are two types of NN interactions. 
The first type takes place in the triangles of atoms 1-2-3 (4-5-6), which are either expanded 
(the case of YMn0 3 ) or contracted (the case of LuMn0 3 ). The second type takes place in 
the bonds 2-1', 2-3', 2-1", 2-3", which are all equivalent. Then, due to the mirror reflec- 
tion x— Y— x, the NN bonds 2-4 and 2-6 between adjacent xy-planes are also equivalent, and 
differ from the bond 2-5'. The same situation holds for the next-NN interactions between 
the planes: there are two equivalent bonds 2-4' and 2-6', which differ from the bond 2-5. 
For the NN interactions, both in and between adjacent xy-planes, there is a clear correla- 
tion between the bondlength and the strength of the exchange coupling. For example, in 
the low-temperature structure of YMn0 3 , the triangle of atoms 1-2-3 (4-5-6) is expanded 
(for two inequivalent NN bonds 2-1 and 2-1' in the xy-pl&ne, the ratio of the bondlengths 
is /2i'/^2i=0.961). Therefore, the AFM interaction J^v is stronger than J21.— The same 
tendency holds for the interplane interactions: for two inequivalent NN bonds 2-5' and 2-4 
(^25'/^24 = U.991), the AFM interaction J251 is stronger than J24. In LuMn0 3 , where the tri- 
angle of atoms 1-2-3 (4-5-6) is compressed, the situation is the opposite: hi'/hi— 1-016 and 
hs'/h^— 1-003. Therefore, the exchange interactions in the bonds 2-1 and 2-4 are stronger 
than in the bonds 2-1' and 2-5'. 
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The behavior of next-NN interactions between t he p lanes obeys quite different rules. 



Since the direct transfer integrals are small (see Ref. [17| for details), these interactions are 
realized as the "super-superexchange" processes via intermediate sites in the pathes 2— >6— >5, 
2— > 1—7-5, etc., which always include one compressed and one expanded bond. Therefore, the 
simple analysis in terms of the bondlengths is no longer applicable. Instead, we have found 
that for all considered compounds (and all considered structures), the AFM interaction in 
the bond 2-5 appears to be weaker than in the bonds 2-4' (and in the equivalent to it bond 
2-6'). Such a behavior has very important consequences: in a noncollinear structure, it is 
more favorable energetically to form the FM coupling in the bond 2-5 in order to maximize 
the AFM coupling in two other next-NN bonds 2-4' and 2-6'. Particularly, it explains why 
the magnetic ground state of YMnC>3 and LuMnOs should be T 3 , T 4 , or r 5 , which are 
characterized by the FM coupling in the bond 2-5, and not T\, r 2 or Tq, where this coupling 
is AFM (see Fig. [2]). In LuMn03, this effect is additionally enhanced by the NN interactions 
between the planes: since AFM interaction in the bond 2-5' is weaker than in two equivalent 
bonds 2-4 and 2-6, it is more favorable energetically to form the FM coupling between the 
sites 2, 5 and 5', where the latter two are connected by the translation. However, in YMnC>3, 
the situation is the opposite and there is a strong competition between NN and next-NN 
interactions between the planes. Particularly, it explains a small energy difference between 
configurations r 3 and T 2 . 

The reliability of the obtained parameters can be checked by calculating #cw- I n the 
case of classical Heisenberg model, the latter is given by the formula 9 CW w £\ J2i/3k B , 
which yields —562 and —650 K for the 10 K structure of YMnC>3 and LuMnC>3, respectively. 
In the quantum case, these values should be additionally multiplied by (1+1/S 1 ). The 
structural changes have some effect mainly on YMnC>3, and, if one uses parameters obtained 
for the 300 K structure, |#cw| decreases by 7% (for comparison, similar change of #cw for 
LuMnC>3 is about 1%). In any case, the obtained values are in a good agreement with 
experimental data.-^^ The calculations of Tn are not straightforward: due to the quasi- 
two-dimensional character of isotropic exchange interactions, Tn will be strongly suppressed 
by thermal fluctuations, as one of the consequences of the Mermin- Wagner theorem.— Of 
course, the molecular-field approximation will overestimate Tn (by factor 4, in comparison 
with the experiment). 

The LMTO calculations yield the following values of NN interactions in the plane xy 
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(in meV): {J 21 ,J 2 i>)= (-12.8,-19.6), (-18.0,-17.0), (-15.8,-16.6), and (-17.0,-17.0) 
for YMn0 3 (10 K), LuMn0 3 (10 K), YMn0 3 (300 K), and LuMn0 3 (300 K), respectively. 
Thus, all the interactions are weaker than in the model analysis. Nevertheless, this seems 
reasonable. First, the NN interactions are generally weaker in the AFM tititi configu- 
ration. This effect was also found in the model calculations, as will become clear below. 
Second, the ratio of AFM to FM contributions to the exchange coupling in manganites 
scales with the value of U as (U— Jk)/(U+3Jh) ~ 1— AJ-r/U.— Thus, larger U, which was 
used in the model (but not in the LMTO calculations), will shift this balance towards the 
AFM coupling. Similar tendency was found for inter-layer interactions: although LSDA, 
supplementing the LMTO calculations, somewhat overestimates FM contributions to the 
exchange interactions, the modulation of these interactions, caused by the lattice distortion, 
again favors the formation of magnetic configurations T 3 or T^. For example, in YMn0 3 
(10 K), the LMTO calculations yield: J 24 = 0.20 meV, J 25 '= 0.04 meV, J 24 '= -0.14 meV, 
and J 2 5= 0.08 meV. Therefore, these calculations confirm that the experimental coupling 
between hexagonal layers is stabilized by the next-NN interactions J25 > ■ The NN 
interactions act in the opposite direction: J 2 ^ < J24. However, their effect is smaller. 

Let us discuss the behavior of the single-ion anisotropy tensor. Due to the mirror reflection 
x— >— x, the tensor t 2 at the site 2 (see Fig. [1]) has the following form: 

^ t xx \ 

f 2 = T yy T yZ 

T Zy T zz J 

where r zy = r yz and T xx +r yy +r zz = 0. Thus, the magnetic moments can either lie along the 
x-axis or be perpendicular to it. In the latter case (and if r yz 7^0) they from a canted magnetic 
structure. The anisotropy tensors at other Mn-sites can be generated by applying the sym- 
metry operations of the space group P6 3 cm. The matrix elements of f 2 can be evaluated in 
the second order of perturbation theory expansion with respect to the SO interactions.— 
Then, near the FM state, we obtain the following sets of independent parameters (in 
meV): (r yy ,r yz ,t zz )= (-0.34,-0.12,0.58), (-0.29,-0.11,0.58), (-0.25,-0.12,0.57), and 
(-0.26,-0.12,0.57) for YMn0 3 (10 K), YMn0 3 (300 K), LuMn0 3 (10 K), and LuMn0 3 
(300 K), respectively. Since r zz >r yy , all structures with large ^-components of the magnetic 
moments are energetically unfavorable. Then, by diagonalizing f 2 , one can find that the 
lowest-energy configuration in LuMn0 3 is the one where the magnetic moment at the site 
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2 is parallel to the x-axis. The next, canted magnetic configuration, is higher in energy by 
about 0.05 meV (for the 10 K structure). This situation is reversed in YMnC>3, where the 
lowest energy corresponds to the canted magnetic configuration. The angle a, formed by the 
magnetic moment and the y-axis, is about 7°. In the next configuration, which is higher in 
energy by about 0.10 meV (for the 10 K structure), the magnetic moment is parallel to the 
x-axis. This energy difference is reduced till 0.01 meV for the 300 K structure. The same 
behavior was found in the LMTO calculations: for YMn03, the lowest energy corresponds 
to the canted magnetic configuration (the canting from the y-axis is about 6°). The next 
configuration, where the magnetic moment is parallel to the x-axis, is higher in energy by 
0.09 eV for the 10 K structure, and this energy difference further decreases for the 300 K 
structure. 

Thus, the change of the ground state from T 3 to T 4 in the direction from YMnC>3 to 
LuMnOs is related to the behavior of the single-ion anisotropy, which in turns correlates with 
the distortion of the 1-2-3 triangles (expansion and contraction, respectively). Moreover, due 
to the 180° rotation around the z-axis, which is required in order to transform the site 2 to 
the site 5 (see Fig. [1]), the matrix element r yz will change sign. Therefore, the canting of 
spins in the planes z=0 and z=c/2 of the 1^ structure will act in the opposite directions, 
and the magnetic moments along the z-axis will cancel each other. 

The single-ion anisotropy will tend to align z-components of the magnetic moments fer- 
romagnetically in each of the xy-pl&ne. However, this effect will compete with the NN AFM 
interactions J21 and J21'. The corresponding analytical expression for the spin canting can 
be obtained by minimizing the energies of single-ion anisotropy and isotropic exchange in- 
teractions: by assuming that all neighboring spins in the xy-plane form the 120°-structure 
(as in the case of the T 2 and T 3 configurations), one can find that 

2t v * 

tan 2a = — , (3) 

T yy _ T zz + 3 J 21 + 6 J 21 , ' V ) 

where the minus-sign corresponds to the situation, which is realized in our HF calculations 
and where e 2 is antiparallel to the y-axis (see Fig. [2]). Then, for the T3 configuration of 
YMnC>3 (10 K), the canting angle a can be estimated (using both model and LMTO param- 
eters of magnetic interactions) as a ~ — T yz ji^J%\ + 6J21') = —0.03°, which is about 7 times 
smaller than the values obtained in self-consistent HF calculations (Table H]). Nevertheless, 
there is an additional contribution to the spin canting, caused by the DM interactions. 
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Parameters of DM interactions between NN sites in the xy-plane are listed in Table HVl 
They were obtained by considering mixed perturbation theory expansion with respect to 



TABLE IV. Parameters of Dzyaloshinskii-Moriya interactions (measured in meV), calculated in 
the ferromagnetic states of YMnOs and LuMn03. The atomic positions are explained in Fig. [TJ 
Calculations have been performed using the experimental parameters of the crystal structure, 
measured at 10 K and 300 K (as denoted in the notations). 

bond YMn0 3 (10 K) LuMn0 3 (10 K) YMn0 3 (300 K) LuMn0 3 (300 K) 

2-1 (0.01,0.01,0.20) (0.04,0.02,0.25) (0.04,0.02,0.17) (0.07,0.04,0.25) 

2-1' (0.03,-0.02,0.21) (0.03,-0.02,0.26) (0.03,-0.01,0.18) (0.02,-0.01,0.26) 

2-1" (0,0.04,0.21) (-0.01,0.04,0.26) (0,0.03,0.18) (0,0.02,0.26) 



the SO interaction and infinitesimal rotations of spin magnetic moments.— In principle, the 
parameters d 2 i' and d2i» are not independent and can be transformed to each other using 
symmetry operations of the space group P6 3 cm. However, it is more convenient to consider 
their contributions independently. Due to the mirror reflection x—>—x, the elements of two 
axial vectors d 23 and d 2 i (see Fig. [T|) obey the following rules: d| 3 = <^ 23 = — ^21 > an d 
d\ z = —d?>i (similar situation holds for other NN interactions). Thus, they will produce a 
finite canting at the site 2 only if the directions of two other magnetic moments e 2 and e 3 
would have an AFM component along x and a FM component along y, i.e.: e\ = — ef and 
e| = e\ . Such a situation is realized in the magnetic configurations T 2 and T 3 (but not in Ti 
and T^). Then, the magnetic moment at the site 2 will experience the additional rotational 
force from the sites 1, 1', and 1": fi_>2 = [d 2 i x ei] + [d 2 i' x ei] + [d 2 i// x ei]. For the magnetic 
configurations T 2 and T 3 , the sites of the type '3' will create the same rotational force: 
f3^2 = However, for the Ti and T 4 configurations, it holds f 3 _> 2 = — fi_> 2 . Therefore, 

these two contribution will cancel each other and there will be no canting of spins. 

These rotational forces should be incorporated in the expression ([3]) for the spin canting, 
which yields a ~ — [r yz + /jl^ 2 )/(3J 2 i + 6J21') = —0.04°. This canting is still smaller than 
a ~ —0.21°, obtained in the HF calculations for the T 3 configuration (Tabled]). Nevertheless, 
it should be noted that all the parameters of the spin Hamiltonian (j2J) were evaluated using 
perturbation theory expansion near the collinear FM state, which is very far from the ground- 
state configuration T 3 . Thus, it is difficult to expect that the perturbation theory, although 
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is very useful for the semi-quantitative analysis, should be able to reproduce all details of 
the solutions of the electronic model ([1]). In fact, some parameters of the spin Hamiltonian 
d2J) appear to be sensitive to the state, in which they are calculated. For example, we 
have also considered the collinear AFM configuration titi'No where the arrows stand for 
the directions of magnetic moments at the sites 1-6. In this case, the DM interactions 
involving the site 2, which is AFM coupled with all NN spins in the xy-pl&ne, become (in 
meV): d 2 i = (0.01,0.03,0.01), d 2V = (0.04,0,0), and d 2 i« = (-0.01,0.05,0.01). Then, 
corresponding rotational force /i_>. 2 wn ^ be about 2 times larger than in the FM state. 
Meanwhile, the parameters of isotropic exchange interactions J 2 i and J%v decrease by about 
15%. These factors will additionally increase a. 

Furthermore, the HF potential for the low-energy model ([1]) is orbitally dependent. In this 
cases, the local force theorem is no longer valid.— Therefore, the total energy change due to 
the SO interaction can be replaced only approximately by the change of the single-particle 
energies of the HF method. For the single-ion anisotropy, the situation was discussed in 
Appendix B of Ref. |23|. Presumably, this is the main reason, explaining the quantitative 
difference between the results of the electronic and spin models. Thus, these are typical 
uncertainties, supplementing the construction and analysis of the spin model (T5]). 

Nevertheless, the local force theorem is valid within LSDA. Therefore, it is interesting 
to estimate the spin canting in the LMTO calculations, which are based on the LSDA 
functional. In this case, all DM interactions become larger. For example, for YMnC>3 
(10 K) we have obtained the following parameters (in meV): d 2 i = (—0.01,0.14,-0.20), 
d 21 , = (-0.16,0.04,-0.12), and d 21 » = (0.06,0.18,-0.26). Then, by combining them with 
corresponding parameters of the single-ion anisotropy r yz = —0.078 meV and isotropic 
exchange interactions J 2 i an d J21', which are listed above, we obtain the canting angle 
a = —0.12°. Thus, it is interesting that LSDA, despite its limitation, provides the best 
starting point for the analysis of the spin canting via the perturbation-theory expansion for 
the spin-orbit interaction, due to validity of the local force theorem. Similar situation was 
found in the orthorhombic LaMnC>3.— 

Thus, although derivation of parameters of spin model (j2j) may differ in details, depending 
on the form of the electronic Hamiltonian, which is used as the starting point, as well as some 
additional approximations, underlying definitions of the model parameters, this analysis 
provides a clear microscopic basis for understanding the main difference between YMnC>3 
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and L11M11O3: why the former tends to form the canted noncollinear magnetic structure T3, 
while the latter forms the planar structure 

C. Magnetic Contribution to Ferroelectric Polarization 

Finally, we would like to comment on the behavior of electronic polarization P||c. 
It was calculated within the Berry-phase formalism, 27 which was adopted for the model 
calculations.- Of course, the ferroelectric activity in YMnC>3 and LuMnOs is primarily 
caused by structural effects. For example, in YMnC>3, the ferroelectric transition occurs 
at about Tq = 880 K,- which is much higher than Tn = 75 K.- This fact was also con- 
firmed by first-principles calculations.— Another appealing evidence is that the ferroelectric 
domains in YMJ1O3 always coincide with the structural ones.- Nevertheless, beside this 
structural deformation, we have found that there is a substantial magnetic contribution to 
P||c. More specifically, all magnetic configurations can be divided in two group. The first 
one includes T\, r 2 , and Tq, where the magnetic moments in the planes z=0 and z=c/2 
can be transformed to each other by the simple rotations. The second groups includes r 3 , 
T 4 , and T 5 , where these rotations should be additionally combined with the time inversion. 
According to our finding, the states in each group are characterized by nearly equal values 
of P||c. However, the transformation of the magnetic state from one group to another 
would cause a finite jump of electronic polarization. Thus, in principle, the value of the 
ferroelectric polarization can be controlled by changing the magnetic state (and vice versa). 
In this sense, more promising candidate is YMnC>3, where the ground state (r 3 ) and the 
first excited state (r 2 ) belong to different groups. The energy difference AE between these 
two configurations is about 0.16 meV (see Table HT1) . Then, the change of the ferroelectric 
polarization, associated with the change of the magnetic state r 3 — KT 2 , can be estimated 
as AP||c= —120 /zC/m 2 . The practical realization of such a switching phenomenon would 
be probably interesting, although it is not immediately clear, which external interaction 
could switch the magnetic state. Formally speaking, the magnetic configuration T 2 could 
be stabilized by the external electric field E||c, which couples to AP and results in the 
additional energy gain —APE. Alternatively, one could exploit the fact that r 2 allows for 
a weak ferromagnetism along z (while T 3 does not) and, therefore, could be also stabilized 
by interaction with the external magnetic field, — AMB, which couples to the net magnetic 
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moment AM (~ — 0.01/iB per Mn-site). However, in order to overcome the total energy 
difference AE, this would require unrealistically large values of E and B, which cannot be 
realized in practice. Therefore, one should explore alternative possibilities. For example, 
from the viewpoint of microscopic interactions, one could use the competition of the NN and 
next-NN interactions between adjacent xy-planes, which in the case of YMnC>3 act in the 
opposite direction (see discussions above). The 1^ configuration is stabilized by the next- 
NN interactions. However, if one could find such macroscopic conditions, which would shift 
this balance towards NN interactions, one could switch the magnetic structure I^— 7T2 and, 
therefore, the ferroelectric polarization. Another possibility is, of course, to exploit mag- 
netism of the rare-earth ions, which can act similar to B, but produces much stronger effect 
on the Mn-sublattice. Such a magnetic phase control was indeed realized experimentally in 
the series of hexagonal manganites with the magnetic rare-earth sublattices.- 1 ^ 



IV. SUMMARY 

Using results of first-principles electronic structure calculations, we have established the 
low-energy model, which is able to capture basic magnetic properties of hexagonal man- 
ganites. This Hubbard-type model describes the behavior of the Mn 3d bands, being sub- 
jected to the lattice deformation and on-site electron-electron interactions. All parameters 
of such model, derived from the first principles calculations for two characteristic manganites 
YMnC>3 and LuMnOs, are summarized in Ref. I7I 

Then, the model was solved in the HF approximation, by considering all possible non- 
collinear magnetic structures with different symmetries. Since the magnetic frustration in 
the hexagonal PQ^cm lattice is lifted by the relativistic SO interaction, the HF approxi- 
mation provides a good starting point for the analysis of the magnetic properties of these 
compounds and successfully reproduce the experimental change of the magnetic ground state 
from PQ 3 cm to PQ^cm in the direction from YMnC>3 to LuMnOs, which was observed in the 
neutron diffraction and nonlinear optical studies. 

In order to clarify the microscopic origin of such a change, we have further transformed 
the electronic model into the spin one and discussed the same trend in terms of differences 
in the behavior of magnetic interactions in these systems. We have found that the main 
reason why YMnC>3 and LuMnOs tend to form different magnetic structure is related to 
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the behavior of the single-ion anisotropy, which couples to the trimerization distortion in 
the hexagonal plane and reflects different directions of this trimerization in YMnC>3 and 
LuMn03 (expansion and construction of the Mn-trimers, respectively). On the other hand, 
the interplane coupling in both compounds is controlled by the next-NN interactions, which 
is less sensitive to the direction of trimerization. The spin canting in the P6gcm structure 
of YMnC>3 is a combined effect of both single-ion anisotropy and Dzyaloshinskii-Moriya 
interactions, which act in the same direction. As the trimerization distortion decreases 
with the temperature, all anisotropic interactions also decrease, thus reviving the magnetic 
frustration and the degeneracy of the magnetic state. 

Finally, using the Berry-phase formalism, we have estimated the magnetic contribution 
to the ferroelectric polarization and discussed how it can be controlled by changing the 
magnetic structure of YMnC>3. 
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